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Abstract 

CP asymmetries in B^ — > D + D~ , D* + D~ , D + D*~ and D* + D*~ decays are inves- 
tigated with the help of the factorization approximation and isospin relations. We find 
that the direct CP violation is governed only by the short-distance penguin mechanism, 
while the indirect CP asymmetries in Bd — > D^D*^ transitions may be modified due 
to the final-state rescattering effect. An updated numerical analysis shows that the 
direct CP asymmetry in vs B® — > D + D~ decays can be as large as 3%. The CP- 
even and CP-odd contributions to the indirect CP asymmetry in B^ vs B® — > D* + D*~ 
decays are found to have the rates 89% and 11%, respectively. Some comments on the 
possibilities to determine the weak phase (5 and to test the factorization hypothesis are 
also given. 
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1 Introduction 



A direct measurement of the CP- violating parameter sin 2/3 in B^ vs P° — > J/ipKs decays, 
where (3 = arg[— (V* b V t d) / (V* b V c d)} is known as an inner angle of the Cabibbo-Kobayashi- 
Maskawa (CKM) unitarity triangle 

v: b v ud + v; b v cd + v*v td = o , (i) 

has recently been reported by the CDF Collaboration [[jj]. The preliminary result sin 2/3 = 
0.79^q;44 (stat + syst) is consistent very well with the standard-model prediction for sin 2/3, 
obtained indirectly from a global analysis of current data on B^-B® mixing, and 

CP violation in K°-K° mixing ||. If the CDF measurement is confirmed, CP violation of 
the magnitude sin 2(3 should also be seen in the decay modes B® vs B® — > D + D~, D* + D~, 
D + D*~ and D* + D*~ , whose branching ratios are all anticipated to be of O(10 -4 ). Indeed the 
channel B^ — > D* + D*~ has been observed by the CLEO Collaboration 0, and the measured 
branching ratio B(D* + D*~) = [6.2^9 (stat) ± 1.0 (syst)] x 1CT 4 is in agreement with the 
standard-model expectation. Further measurements of neutral and charged B decays into 
£)WdW states will soon be available in the first-round experiments of KEK and SLAC B- 
meson factories as well as at other high-luminosity hadron machines (see, e.g., Ref. |4j] for a 
review with extensive references). 

In the literature some special attention has been paid to B — > £)(*)/)(*) transitions and 
CP violation. For example, the CP properties of B d — > _D(*)+.D(*)~ decays were analyzed in 
the heavy quark limit in Ref. ||; the isospin relations and penguin effects in B — > £)(*).£)(*) 
decays were explored in Ref. ||; the possibility of extracting the weak phase (3 and testing 
the factorization hypothesis in B^ vs B^ decays into the non-CP eigenstates D ± D* T were 
investigated in Ref. 0; and the angular analysis of B d — > D* + D*~ decays to determine 
CP-even and CP-odd amplitudes were presented in Ref. [^]. In addition to those works, 
numerical estimates of branching ratios and CP asymmetries in B — > D^D^ decays have 
been given in Ref. 0, in which neither electroweak penguin contributions nor final-state 
rescattering effects were taken into account. 

The present paper, different in several aspects from those previous studies, aims at ana- 
lyzing final-state rescattering effects on direct and indirect CP asymmetries in B — > D^*'D^*> 
decays. We calculate the 1 = 1 and 1 = isospin amplitudes of these processes by using 
the factorization approximation and the effective weak Hamiltonian, and account for long- 
distance interactions at the hadron level by introducing elastic rescattering phases for two 
isospin channels of the final-state mesons. In this approach we find that direct CP asym- 
metries in both charged and neutral B decay modes are governed only by the short- distance 
penguin mechanism, but indirect CP asymmetries in B d — > D ± D* T transitions may be 
modified due to the final-state rescattering effect. An updated numerical analysis of direct 
CP violation in B — > DD, D*D, DD* and D*D* decays is made without neglect of the 
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electroweak penguin effects. We obtain the asymmetry as large as 3% in P+ — > D + D° vs 
P~ — > D~D° or 5° vs P^ — > D + D~ decays. In the absence of angular analysis we find 
that the indirect CP asymmetry in — ► D* + D*~ decays is diluted by a factor 0.89, i.e., 
11% of the asymmetry arising from the P-wave (CP-odd) contribution. We also give some 
comments on the possibilities to determine the weak phase j3 and to test the factorization 
hypothesis in the presence of final-state interactions. 



2 Isospin amplitudes 



The effective weak Hamiltonian responsible for B — > p(*)p(*) decays can explicitly be written 
as @ 

/ 2 10 

v gb v* d £ Ci Qi + y, aQi 



K 



cff 



v/2V 
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Wilson coefficients, and 
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u, c) are the CKM matrix elements, q (for z 
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10) are the 



Q2 = (<*7) V -a(?&)v-a ' 
Qs = (*)v-a(cc)v-a » 

C/4 = (^a&/3)v-A(C/3C a ) V --A 5 

Q 5 = (*)v-a(cc)v+a » 

Qe = (d a b/3) v _ A (ci3C a ) v+A , (3) 

as well as Qi = Q5, Qg = Qq, Qq = Q3 and C/io = Qi- Here Q3, • • Q§ denote the QCD- 
induced penguin operators, and Q 7 , • • •, Q w stand for the electroweak penguin operators. 

It is clear that the AB = +1 and AB = — 1 parts of H, e s have the isospin structures 
|l/2, +1/2) and |l/2, -1/2), respectively. They govern the transitions P+ -> pM+P)W°, 
P° — > p(*)°Z)(*) and their CP-conjugate processes. The final state of 



£0 £)(*)+£)(*)- 



each decay mode can be in either I = 1 or / = isospin configuration. For simplicity we 
denote the amplitudes of six relevant transitions by use of the electric charges of their final- 
and initial-state mesons, i.e., A +0 , A + ~ , A 00 (for P + and P° decays) and A~°, A + ~ , A 00 (for 
P~ and P^ decays). These amplitudes can be expressed in terms of the 1 = 1 and 1 = 
isospin amplitudes, which include both weak and strong phases. For example R, 



A+ o 

A + - 



A 1 
1 



(A, + A Q ) 



(A, - A 







(4) 



2 As for B — > D*D* decays, the isospin relations hold separately for the transition amplitudes with helicity 
A = -1, or +1. 
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and the relations between (A °, A + , A 00 ) and (A\, Aq) hold in the same form. In the 
complex plane two sets of isospin relations form two triangles: A +0 = A + ~ + A 00 and 
A~° = A+- + A 00 . 

To calculate the magnitudes of / = 1 and 1 = isospin amplitudes, we make use of the 
effective Hamiltonian 7i e g and the factorization approximation. We neglect the contribu- 
tions of the annihilation-type channels, which are expected to have significant form-factor 



suppression It should be noted that in this approach the Wilson coefficients and the 
relevant hadronic matrix elements of four-quark operators need be evaluated in the same 
renormalization scheme and at the same energy scale. Following the procedure described 



in Ref. |12[ one can obtain the scale- and renormalization-scheme-independent transition 
amplitudes consisting of the CKM factors, the effective Wilson coefficients, the penguin loop- 
integral functions and the factorized hadronic matrix elements. Under isospin symmetry, we 
are only left with two different hadronic matrix elements: 

Z = ( J D W+ |(crf) V -A|0)( J D W -|(&c) v _ A |5°) 

= ( J DW + |(crf) v _ A |0)( J DW°|(6c) v _ A |5 u + ), 

Z = ( J D W -|(rfc) v _ A |0)( J DW + |(c6) v _ A |5 rf °) 

= ( J D W -|(rfc) v _ A |0)( J DW°|(c6) v _ A |5-). (5) 

Note that \Z\ = \Z\ holds for the final states with two pseudoscalar mesons or those with 
one pseudoscalar and one vector mesons. Only for the final states with two vector mesons 
\Z\ and \Z\ are different, as the P-wave contributions to Z and Z have the opposite signs 
(see section 4 for the detail). Furthermore, we account for final-state interactions at the 
hadron level by introducing the elastic rescattering phases S\ and Sq for 1 = 1 and 1 = 



isospin channels (a similar treatment can be found, e.g., in Refs. []13 , |i"4|). We then arrive 
at the factorized isospin amplitudes as follows: 

a x = ^( Vud v: b s u + v cd v; b s c )Ze iS \ 

A = % (V ud V* b S u + V cd V* b S c ) Ze i5 ° , (6) 

in which S u and S c are composed of the effective Wilson coefficients and the penguin loop- 
integral functions (see section 3). The expressions of A\ and A can be obtained respectively 
from those of Ai and A in Eq. (6) through the replacements Z =^> Z and V qd V* b V* d V gb 
(for q = u and c). Note that all parameters in the isospin amplitudes, except the CKM 
factors, are dependent upon the specific final states of B decays. 

One can see that \A$\ = \Ai\ and \Aq\ = \A\\ hold in the context of our simple factoriza- 
tion scheme. This implies that the B d — > Z)(*) z)(*) transitions would be forbidden, if there 
were no final-state rescattering effects (i.e., if 5 = Si). Substituting Eq. (6) into Eq. (4), 
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one obtains 



A + - = ^{V u ^: b S u + V c ^: b S c )ZcoJ-^e^ + ^\ 

Similarly A + ~ and A 00 can be read off from A + ~ and A 00 through the replacements Z =>- Z 
and V q dV* b =>- V* d V q b (for g = -u and c). It is easy to find 



\A + - 


2 _|_ ^00 2 _ 


\A +0 \ 2 


\A + - 


2 _|_ ^00 2 _ 


\A- Q \ 2 



i.e., the two isospin triangles are right-angled triangles. Whether the relations in Eq. (8) 
are practically valid or not can be checked, once the experimental data on branching ratios 
of B — > £)(*)£)(*) decays are available. 

If \A + ~\ = \A +0 \ held, \A 00 \ = would result within the factorization approach described 
above. Namely, observation of the (approximate) equality between the decay rates of B d — > 
£)(*)+£)(*)- anc i B+ -> D(*)+5W° would imply that the decay modes 52 -> D^°D^° were 
forbidden or strongly suppressed. This conclusion is in general not true, however. Without 
any special assumption or approximation, we denote Aq/Ai = ze 10 , \A 00 /A +0 \ 2 = R and 
obtain consequences of the equality \A + ~\ = \A +0 \ as follows: 

z = V3 + cos 2 6 — cos 6 , 

R = 1 + cos 2 9 - cos #V3 + cos 2 6 . (9) 

The behaviors of z and R changing with 9 is illustrated in Fig. I. It is clear that in general 
I A 00 1 = (i.e., R — 0) is not necessary to hold even if \A + ~\ = \A +0 \ holds. Therefore 
the detection of B d — > /}(*)°/)(*) transitions is very useful in experiments, in order to 
demonstrate whether final-state rescattering effects are significant and to test whether the 
factorization approximation works well. 



3 Direct CP asymmetries 

We proceed with the factorization scheme to calculate direct CP asymmetries in the decay 
modes under discussion. As for the final states with two vector mesons, we sum over their 
polarizations and arrive at \Z\ 2 = \Z\ 2 , a relationship which apparently holds for other types 
of final states. With the help of Eqs. (6) and (7) it is easy to show that the decay rate 
asymmetry between B+ — ► Z)(*)+Z)(*) and B~ — ► £)(*) _ £)(*) decays is identical to that 
between B Q d and B° d -> £>«+£>(*)- or L>(*) Z)(*) decays. All these CP asymmetries are 
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Figure 1: Behaviors of z and R changing with 9 under the condition \A + | = |^4 +0 |. 
independent of the rescattering phases and the hadronic matrix elements Q 

A = 
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(10) 



L4 00 | 2 + |A 00 | 2 

2rsin 7im(£ u £* 



r 2 |C| 2 +lCc| 2 -2rcos 7 Re(C tt C*) ' 

where r and 7 are defined by re 17 = — (V u dV* b ) / (V c dV* b )- The phase 7 corresponds to another 
inner angle of the unitarity triangle defined in Eq. (1). Eq. (10) indicates that direct CP 
violation arises only from final-state interactions of the quark level (through the penguin 
mechanism) in B — > D^*'D^*> decays. This result, as a straightforward consequence of the 
factorization approximation, can directly be confronted with the upcoming experiments at 
S-meson factories. 

Let us evaluate the direct CP asymmetries A for different final states. As mentioned 
above, S u and S c in Eq. (10) depend on the effective Wilson coefficients q and the penguin 
loop-integral functions F q . The latter can be given, for a momentum-squared transfer k 2 at 



3 For time-integrated Bd decays the direct CP asymmetries are diluted by a well-known factor 1/(1 +x d ), 
where Xd ~ 0.7 is the B^-B® mixing parameter. In this paper we do not take such mixing effects into 
account. 
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the 0(rrib) scale, as follows ITS 



4 dx x(l — x) In 



m — k x(l 

mi 



x) 



'IV. 



The absorptive part of F q , which is a necessary condition for direct CP violation, emerges 
if k 2 > 4m 2 ,. The concrete expressions of S u and S c are found to be 
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where 



Cl — "g - + C 4 + + , 

C 2 = j + cz + Ci, 

= y + c 6 + — + c 8 , 

C 4 = c 2 a s + (c! + «e 



(13) 



In these equations q (for i = 1, • • •, 10) are the renormalization-scheme- independent Wilson 
coefficients, a s and a e stand respectively for the strong and electroweak coupling constants, 
and £ is a factorization parameter arising from the transformation of (V-A)(V+A) cur- 
rents into (V-A)(V-A) ones for the penguin operators Q5, • • •, Q$- Note that £ depends on 



properties of the final-state mesons |TB 

£ = 



2m 2 D 



(DD) 



(m c + m d )(m b -m c ) 

(D*D) 

D 



2m 2 , 



(m c + m d ){m b + m c 




(DD*) 
(D*D*) 



(14) 



where the order of two D^ mesons corresponds to that in the factorized hadronic matrix 
element Z or Z, as given in Eq. (5). 

With the help of Eqs. (11) - (14) we are able to calculate the CP asymmetries A 
numerically. Note that \S U \ \S C \, as the former consists only of the penguin contribution 
and the latter is dominated by the much larger tree-level contribution. This, together with 
|r| < 1, allows an instructive analytical approximation of A: 



A ~ 2r sin 71m ( ^ 



(15) 



For illustration, we typically choose m u = 5 MeV, m c = 1.35 GeV, m b = 5 GeV and m t = 174 
GeV. The strong coupling constant is taken as a s = 0.21 at the 0(m b ) scale. Values of 



7 



1 1 1 1 r 

A{DD) 



A(DD*) x 10 



A(D*D) = A(D*D*) 



j I 

0.2 0.4 0.6 0.8 1 

k 2 /ml 

Figure 2: Direct CP asymmetries of B — > D^D^*' in the factorization approximation. 

the effective coefficients c* read |T7|1 : c\ = —0.313, C2 = 1.150, C3 = 0.017, C4 = —0.037, 
c 5 = 0.010, c 6 = -0.046, c 7 = -0.001a e , c 8 = 0.049a e , c 9 = -1.321a e and c 10 = 0.267a e 
with a e = 1/128. The CKM factors are taken to be r = 0.38 and 7 = 60°, consistent with 
the lastest data on quark mixing and CP violation ||. The unknown penguin momentum 
transfer k 2 is treated as a free parameter changing from O.Olm^ to m\. Our numerical results 
are shown in Fig. 2. Some discussions are in order. 

1. All CP asymmetries have the same sign and undergo a change of magnitude at k 2 = 
4m 2 ~ 0.3?7i|. The asymmetry A(DD) is most sensitive to the uncertain penguin 
momentum transfer k 2 , but its magnitude increases only about 0.3% from k 2 = O.Olm^ 
to k 2 = ml. It is found that the strong (gluonic) penguin effect is dominant over the 
electroweak penguin effect, thus the latter is safely negligible. 



0.03 



0.025 



0.02 



0.015 
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2. The CP asymmetry A(DD) can be as large as 3%, while A(DD*) is only about 
2 x 10~ 3 . The smallness of the latter comes from the cancellation effect, induced by 
the factor (1 + £) with £ ~ —0.8, in S u and S* c . In our factorization approximation, 
the asymmetries A(D*D) and A(D*D*) are identical and of magnitude 1%. 

3. Observation of the CP asymmetries A{DD) and A{D*D*) to three standard deviations 
needs about 10 8 B^ events, if the composite detection efficiency is at the 10% level. 
More events are required to measure the same CP asymmetries in Bd decays, due to 
the cost for flavor tagging. 

It is therefore worth while to search for such direct CP- violating signals in the first-round 
experiments of P-meson factories. 



4 Indirect CP violation 



Although the final-state rescattering phases have no effect on direct CP asymmetries A in 
our factorization scheme, they are possible to influence the indirect CP violation arising from 
the interference between direct Bd transition and B^-B® mixing in the decay modes under 
consideration. The characteristic measurable of this source of CP violation is in general a 
difference between two rephasing-invariant quantities defined as |Tj| 



A(/) 
*(/) 



Im 


£ 
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Im 


P 




Q 



A(B° d - /) 



A(B° d 
A(B° d 



f) 
I) 



A(B° d - /) 



(16) 



where q/p = (V t lV t d) / (VtbV t * d ) denotes the weak phase of Bj-B® mixing [], and / is the CP- 
conjugate state of /. If / is a CP eigenstate (i.e., |/) = CP\f) = ±|/)) and the decay is 
dominated by the tree- level channel, then A(/) = — A(/) is a good approximation. In general 
only the difference A(/) — A(/), which will vanish if all the CKM factors are real, measures 
the CP asymmetry. Note that the CP-even and CP-odd components of / = D* + D*~ state 
or / = Jj*°l)* state may cause some dilution in the measurables A(/) and A(/). A proper 
treatment of indirect CP violation in such modes is to make use of the angular analysis ||. 
Alternatively one may evaluate the P-wave contribution to A(P* + P*~) and A(P*°P*°) by 
use of the factorization approximation and the heavy quark symmetry ||20|1 . 



As penguin contributions to the transition amplitudes of B — > D^D^*' decays have been 
estimated to be at the percent level, we expect that their effects on A(/) and A(/) are 
unimportant and negligible. 



4 The CP violation induced solely by B^-B^ mixing (i.e., \q/p\ ^ 1) is expected to be negligibly small (of 
order 10~ 3 or smaller |Is|]) in the standard model. 
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It is obvious, as shown in Eq. (7), that for and decays into the CP eigenstates 
D + D~ and D°D° the amplitude ratios A^ / A + ~ and A 00 /A 00 are independent of the rescat- 
tering effects. Neglecting small penguin contributions to S u and S c (i.e., taking S u — and 
S c = c 2 + ci/3), we arrive at 

A{D + D-) = A{D°D°) = + sin 2/5, 

A(D+Zr) = A(D°D°) = -sin 2/3, (17) 

where /3 is just the inner angle of the unitarity triangle defined in Eq. (1). Therefore the 
measurement of indirect CP asymmetries in Bd — > D + D~ and B^ — > D°D° decays may 
serve as a cross-check of sin 2/3 extracted from the CP asymmetry in Bd — > J/ipKs decays. 

The channels B d -> D* + D-, D + D*- and P d -> P* /) , P°P*°, whose final states are 
non-CP eigenstates, are also useful for extraction of the weak angle (3. Since the pseudoscalar 
and vector mesons from P° and those from Bd have different quark-diagram configurations, 
the hadronic matrix elements and final-state rescattering phases in these two processes should 
in general be different 0. As a result, 

A(D* + D-) = CP+- sin(5 + 2(3) , 

A(D + D*-) = CP+- sin(5 - 2/3) ; (18) 



and 



A(P*°P°) = CPoosin(5 + 2/?) , 

A(P°P*°) = CPoosin(5-2/3) , (19) 



where 
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2 

In obtaining these results we have neglected the small penguin effects. The decay constants 
and form-factors in the expression of (, coming from decomposition of the hadronic matrix 
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elements Z d q* and Z D *£, given in Eq. (5), are self-explanatory. Note that P + _ = 1 and 
5 = 5®^* — 5®*® hold, if one takes the limit 5{ = 5$ (for each final state /), in which the 
decay modes Bd — > D^*'°D^*'° become forbidden. In the presence of rescattering effects, i.e., 

P H 7^ 1, the extraction of (3 from A(D* + D~) and A(D + D*~) seems difficult. However, 

it is possible to determine the isospin phase difference S{ — 5$ from the triangle relation in 
Eq. (4), if the relevant rates of three (one charged B and two neutral B) decay modes are 
measured in experiments. The observation of Bd — > D*°D° and D°D*° transitions turns out 
to be crucial: (a) if their branching ratios in comparison with those of Bd — > D* + D~ and 
D + D*~ decays are too small to be detected, then the final-state rescattering effects should 
be negligible and the naive factorization approach with P + _ = 1 might work well; (b) if 
their branching ratios are more or less comparable with those of Bd — > D* + D~ and D + D*~ 
decays, then a quantitative isospin analysis should be available, allowing us to extract the 
isospin phase differences and determine the magnitudes of P+_ and Poo- in both cases, ( can 
experimentally be determined and the result can be confronted with the theoretical value of 
( calculated by inputting relevant decay constants and form-factors. 

For Bd — > D* + D*~ and Bd — > D*°D*° decay modes the indirect CP asymmetries need a 
more careful analysis. Note that the transition amplitude of B® — > D* + D*~ (or is 
a sum of three different components, i.e., the S-, D- and P-wave amplitudes Without 



loss of generality the hadronic matrix elements Z and Z for Bd — > D* + D* can be written 

as 

Z = fi(e+-e_) + -\- (p -e+)(p -e_) 

m 



D* 
m D* 



a (e + -e_) + — (p -e+)(p 

m D* 



- i-5- (e a We+ n e«p^p nX ) , (22) 



m 



D" 

where e± denotes the polarization of D* ± meson, p and p± stand respectively for the mo- 
menta of Bd and -D* ± mesons, and (a, b, c) are real scalars without the penguin effects. In 
terms of the decay constants and form factors, a, b and c read explicitly as 

a = m D *f D *(m B + m D *) A^ D * (m 2 D «) , 



-2m 3 D .f D . 



c = -2m D *f D * ■ . (23) 

m B + m D * 

In the absence of angular analysis one may first calculate the ratio ZjZ by summing over 
the polarizations of two final-state vector mesons [21], and then calculate the CP- violating 
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quantities A(D* + D* ) and A(D* + D* ) in the neglect of small penguin effects. We finally 
arrive at 



A(D* + D*-) = A(D*°D*°) = + sin2/? 

A(D* + D*-) = A(D*°D*°) = -sin2/3^^, (24) 



1 + X 

i-x 

1 + X ' 



where 



X = . (25) 

(2 + x 2 )a 2 + (x 2 - f ) 2 b 2 + 2x(x 2 - 1)56 

with x = (m 2 B — 2m 2 D *) / (2m 2 D ,) = 2.45. Clearly the dilution parameter x results from the 
P-wave contribution to the overall decay amplitudes. If we adopt the simple monopole model 
for relevant form factors [g|], it turns out that V BD * (m 2 D ,) = 0.784, Af D *(m 2 D .) = 0.715 



and A^ * (m 2 D *) = 0.753. Accordingly 6/5 = —0.160 and c/a = —0.167. The relationship 
6/5 ~ c/a is indeed guaranteed by the heavy quark symmetry, which makes the form factors 
appearing in Eq. (23) related to one another. In this symmetry limit we obtain [2C] 



~=~- = r , (26) 

a a m B (m B + 2m D ,) 

amounting to —0.164. Then we get (1 — %)/(l + %) w 0.89, a value deviating only about 11% 
from unity. Note that this dilution factor can also be determined from measuring the ratio 
A(D* + D*~) / A(D + D~) . From this estimation we find that the P-wave dilution effect is not 
very significant, therefore extracting the CP-violating parameter sin 2/3 from Bd — > D* + D*~ 
decays remains possible even if a delicate angular analysis is not made. 



5 Summary 

We have analyzed direct and indirect CP asymmetries in P° vs P° — > D + D~, D* + D~, 
D + D*~ and D* + D*~ decays. The isospin amplitudes of these transitions are calculated with 
the help of the effective weak Hamiltonian and the factorization approximation, and the 
long-distance interactions at the hadron level are taken in to account by introducing elastic 
rescattering phases for two isospin channels of the final-state mesons. We have shown that in 
this factorization approach the direct CP violation is irrelevant to the final-state rescattering 
effects, i.e., it is governed only by the short-distance penguin mechanism. The magnitude 
of direct CP violation is estimated to be 3% in Bd — > D + D~ decay modes. The same 
amount of CP violation can manifest itself in the charged B u decays into D + D° and D~D° 
states, which are easier to be measured at P-meson factories. We have demonstrated that the 
penguin effects on indirect CP asymmetries in Bd — > p(*)+p(*) _ decays are insignificant and 
even negligible. While the long-distance rescattering has no effect on indirect CP violation in 
Bd — > D + D~ and D* + D*~ transitions, it may affect that in Bd — > D ± D* ZfZ modes, whose final 
states are non-CP eigenstates. We have calculated the P-wave contribution to the indirect 
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CP asymmetry in B d — > D* + D*~ decays. The corresponding dilution effect is found to be 
insignificant, therefore observation of large CP violation remains under expectation even 
without the delicate angular analysis. 

It is certainly necessary to test the validity of our factorization hypothesis, on which 
most of the afore-mentioned results depend. To do so a measurement of B d — > £)W°Z)(*) 
transitions will be particularly helpful. On the one hand, if the branching ratios of these 
decay modes are too small comparied with those of B d — > £)(*)+£)(*)- transitions, then the 
final-state rescattering effects should be negligible and the naive factorization approximation 
might work well. On the other hand, if the branching ratios of B d decays into D^°D^° and 
£)(*)+_£)(*)- states are found to be more or less comparable in magnitude, then a quantitative 
isospin analysis should become available, allowing us to extract the isospin phase differences 
and control the final-state rescattering effects. In any case much can be learnt about the 
factorization hypothesis and its applicability in B decays into two heavy charmed mesons. 

In conclusion, the observation of direct and indirect CP asymmetries in B d — > 
decays is promising at 5-meson factories. They are expected to provide us some valuable 
information about the weak phase j3 as well as the penguin and rescattering effects in non- 
leptonic B decays. 
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